Spin-dependent electrical transport in ion-beam sputter deposited Fe-Cr multilayers 
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The temperature dependence of the electrical resistivity 
and magnetoresistance of Xe-ion beam sputtered Fe-Cr multi- 
layers has been investigated. The electrical resistivity between 
5 and 300 K in the fully ferromagnetic state, obtained by ap- 
plying a field beyond the saturation field (Hsat) necessary for 
the antiferromagnetic(AF)-ferromagnetic(FM) field-induced 
transition, shows evidence of spin-disorder resistivity as in 
crystalline Fe and an s-d scattering contribution (as in 3d met- 
als and alloys). The sublattice magnetization m{T) in these 
multilayers has been calculated in terms of the planar and 
interlayer exchange energies. The additional spin-dependent 
scattering Ap(r) = p{T,H = 0)af - p(T,H = H,^t)FM in 
the AF state over a wide range of temperature is found to be 
proportional to the sublattice magnetization, both Ap(r) and 
m(T) reducing along with the antiferromagnetic fraction. At 
intermediate fields, the spin-dependent part of the electrical 
resistivity {ps{T)) fits well to the power law ps(T) = b — cT" 
where c is a constant and b and a are functions of H. At 
low fields a ~ 2 and the intercept b decreases with H much 
the same way as the decreeise of Ap(T) with T. A phase dia- 
gram (T vs Hsut) is obtained for the field-induced AF to FM 
transition. Comparisons are made between the present in- 
vestigation and similar studies using dc magnetron sputtered 
and molecular beam epitaxy (MBE) grown Fe-Cr multilayers. 
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I. INTRODUCTION 

As one of the very few lattice matched transition 
metal pairs one of which is ferromagnetic, Fe-Cr mul- 
tilayers offer excellent opportunities for investigating the 
exchange coupling of Fe layers through an antiferromag- 
netic Cr spacer layer giving rise to the so-called Giant- 
Magnetoresistance (GMR). Applications as magnetic- 
field sensors, especially in reading information, sensing 
position/speed of moving parts, etc. have triggered in- 
tense research activity in these multilayers. GMR sensors 
are not only very sensitive but they can be made very 
small in size. For practical purposes not only does one 



need large GMR but also small saturation fields. 

The aim of the present work is to study the tempera- 
ture dependence of electrical resistivity and magnetore- 
sistance in GMR multilayer stacks prepared by ion-beam 
sputter deposition. Typical multilayers reported here 
comprise 30 repeat layers of [Fe(20A)/Cr(10A)] that have 
been deposited by ion-beam sputter deposition onto Si 
substrates with Xenon ions at 900 V and a beam current 
of 20 mA. 

GMR in multilayers can be understood in terms of 
some simple ideas as follows. In zero magnetic field the 
ferromagnetic Fe layers are coupled antiferromagnetically 
through the Cr spacer layer giving rise to a high electrical 
resistance. This antiferromagnetic coupling is ascribed to 
the indirect exchange interaction between the Fe layers 
through the oscillatory RKKY interaction mediated by 
the conduction electrons. The above antiferromagnetic 
coupling between Fe layers was established by means of 
light scattering from spin waves. 0] As the external field 
increases, the spins in different Fe layers align in the di- 
rection of the field producing a completely ferromagnetic 
alignment beyond a saturation field -ffsat, reducing the 
resistance. Thus we have a negative magnetoresistance 
(MR). 

Magnetoresistance is defined (in text-book fashion) by 



MR = 



p(g,T)-p(0,r) 
P(0,T) 



X 100% 



(1) 



It is found in dc magnetron sputtered Fe-Cr superlattices 
Q that GMR oscillates as a function of Cr spacer thick- 
ness with 3 gradually decreasing peaks centered at 11, 27 
and 42 A Cr for Fe thickness of 32 A. Also the first an- 
tiferromagnetic region occurs between 6 and 11 A Cr for 
Fe thickness lying between 15 and 40 A. We have there- 
fore chosen the Cr thickness around 10 A corresponding 
to the strongest peak of the antiferromagnetic coupling 
between the Fe layers (and hence the highest GMR). 

The basic qualitative features of GMR can be under- 
stood even in terms of bulk scattering only if the mean 
free path (mfp) of electrons within the layers is much 
larger than the layer thickness. If the mfp of the elec- 
trons is larger than the Cr spacer thickness, the electrons 
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can feel the relative orientation of the magnetization of 
the successive layers. However, this interplay between 
the successive magnetic layers disappears and the GMR 
vanishes if the mfp is less than the Cr layer thickness. 

GMR is attributed to the spin-dependent conduc- 
tion properties of ferromagnetic metals. In a ferromag- 
netic metal or alloy the electrical conduction takes place 
through independent channels by spin-up (called major- 
ity) and spin-down (called minority) electrons. This is 
the two-current model of Fert and Campbell |^ whose 
physical basis is the dominance of the spin-conserving 
scattering and the weakness of the spin-flip collision, at 
least at low temperatures. In this picture all electrons 
of a given spin (up or down) with s or d or hybridized 
character are grouped together to form majority (up) or 
minority (down) bands. If one takes into account the 
details of the band structure of Fe (a weak ferromag- 
net) and the simple Drude conductivity formula for each 
band, it is easily shown that the majority band has a 
much higher conductivity than the minority band. As a 
result, in the ferromagnetic alignment brought about by 
the saturation field (-ffsat)j there is hardly any scattering 
for the majority band electrons since they remain major- 
ity in all the Fe layers. On the other hand, the minority 
band electrons get scattered within every Fe layer. Hence 
there is a short-circuiting effect, so to say, and the resis- 
tance (pfm) drops in the ferromagnetic alignment. How- 
ever, in the antiferromagnetic configuration (zero-field) 
both the majority and minority band electrons are scat- 
tered in successive layers, and there is no short-circuiting 
effect. Therefore the resistance (paf) remains relatively 
high. A very simple calculation based on the two-current 
model shows that the GMR at low temperatures is given 
by (considering only bulk scattering) 



GMR: 



piH = i/sat) - p(0) Pfm - PAF 



P(0) 

Pl/Pl - 1 
Pl/Pl + 1 



PAF 



(2) 



where pi and p| are the resistivities of the minority and 
majority carriers, respectively. It turns out that interface 
scattering from imperfect interfaces, defects and impuri- 
ties in Fe-Cr multilayers is also spin-dependent, and gives 
rise to the GMR. As a matter of fact, it is the imbalance 
between the resistivities of the two bands which is re- 
sponsible for the GMR for bulk, interface and spin-flip 
scattering. The subject of GMR has been reviewed very 
well in a recent book. |^ 

Considerable work has been reported on electrical 
transport and magnetic properties of Fe-Cr multilayers 
prepared by sputtering or Molecular Beam Epitaxy 
(MBE). |^-||] The temperature dependence of the elec- 
trical resistivity in [Fe(30A)/Cr(10 - 50A)]xlO multilay- 
ers, prepared by sputtering and MBE, was interpreted by 
Almeida et al. ^ in terms of phonon-assisted s-s and s-d 
scattering in the temperature range of 15-300 K and in a 
saturation magnetic field of 7.5 kOe. 



In the antiferromagnetic Fe-Cr superlattice, made by 
dc magnetron sputtering, the temperature-dependent 
magnetoresistance (defined as pm{T) = p{T,H = 0) — 
p{T,H = Hsat) ) was found by Mattson et al. ||] to fol- 
low the equation for T < 100 K 



PuiT) ^ Pm{T = 0) + ApuiT) 
= Pm{T = 0) - aT^ , 



(3) 



where a is a constant of proportionality. This be- 
haviour was explained in terms of thermal excitation of 
magnons whose occupation number (n) oc (at low 
temperatures) for anisotropic materials H and assuming 
Apm(T) n. 

In [Fc(12A)/Cr(12A)]xlO multilayers, prepared by 
MBE on MgO(lOO) substrates, several interesting obser- 
vations were made by Aliev et al. Among them are: 

(a) The isothermal magnetoresistance, defined as [p{0) — 
p{H)], is proportional to H in the parallel (magnetic field 
in the plane of the multilayers) orientation and to in 
the perpendicular case. 

(b) A T - Hsat phase diagram was obtained, which clearly 
indicated the transition between AF and FM states. 

(c) As opposed to the work of Mattson et al. (which 
looked at the difference between the resistivity in the 
ideally antiferromagnetic (H = 0) and ferromagnetic 
{H = Hsut) alignments) here ^ the spin-dependent part 
of the electrical resistivity, defined as Ps{T) = p{T, H) — 
p{T,H > Hsat) for fields H > (where one has both 
ferromagnetic and antiferromagnetic fractions), is found 
to vary in a wide range of temperature below 100 K as 



Ps{T) 



PsiT^O) 



Ap.(T) 



(4) 



where b = ps(T = 0) and the temperature exponent a 
are functions of the magnetic field H , and c is a constant 
of proportionality. The constant a was found Q to be 
~ 1.7 for ~ 0, ~ 2.0 for H < O.bH^^t, and ~ 1 for 
H ~ Hsut ■ This is in contrast to the value of a = 2 for 
H ~ (purely antiferromagnetic) in the work of Mattson 
et al. II (Eq. (3) above). 

(d) Aps{T) was found to vary linearly with temperature 
from 20 mK to about 1.5 K "which could be due to elec- 
tron scattering on critical thermal spin fluctuations" . j|] 
All the work summarized above is mostly on MBE 
grown Fe-Cr multilayers and those made by dc mag- 
netron sputtering. Although the present work is on Fe- 
Cr multilayers prepared by ion-beam sputter deposition, 
we do not believe that the underlying physics depends 
in any significant way on differences in these deposition 
techniques. 



II. EXPERIMENTAL 

Fe-Cr multilayers were prepared by ion-beam- 
sputtering technique and characterized by Transmission 
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Electron Microscopy (TEM), Atomic Force Microscopy 
(AFM), Auger Electron Spectroscopy (AES), X-ray Pho- 
toelectron Spectroscopy (XPS), as well as resistivity, 
magnetic hysteresis loop and magnetotransport measure- 
ments. The film deposition procedure and film proper- 
ties, including chemical composition, surface morphology, 
resistivity, saturation magnetization, coercive field, and 
magnetoresistance ratio, are discussed in detail in Ref. 
[10]. Values of GMR ratios of the films are comparable to 
values measured for polycrystalline Fe-Cr films deposited 
by the more conventional rf sputtering technique. We 
have deposited the following Fe-Cr multilayer combina- 
tions : 

Si/Cr(50A)/[Fe(20A)/Cr(t(A))]x30/Cr(50-t (A)), where 
t was varied from 8 to 14A; this range surrounds the first 
antiferromagnetic maximum in the Fe-Cr multilayer sys- 
tem. The deposition rates varied from 5 to 30 A/min 
depending on the primary ion beam energy, the type of 
ions, and the target material. The films were deposited at 
room temerature. The effects of variations of the primary 
ion beam energy and the type of ions on GMR values were 
examined; the investigated primary ion energy range was 
700-1200 eV for Ar ions and 900-1200 eV for Xe ions. It 
was demonstrated that the GMR ratio is greater for films 
deposited using Xe ions than for films deposited using Ar 
ions, and that for both types of ions the GMR ratio in- 
creases as the primary ion beam energy decreases. In this 
investigation we report the work on Fe-Cr multilayers of 
typical structure [Fe(20A)/Cr(10A)] x30 layers grown on 
Si substrates using Xe ion at 900 V and a beam current 
of 20 mA. 

The temperature dependence of the resistance between 
5 and 300 K of the Fe-Cr multilayers was measured in 
zero as well as in some applied magnetic fields using the 
standard four-probe dc technique and a magnetic field 
(0-5. 5T) provided by a Quantum Design SQUID mag- 
netometer (MPMS). Both the transport current and the 
applied field were in the plane of the film with the cur- 
rent parallel to the field. For measurements in magnetic 
fields perpendicular to the film plane we used a Quantum 
Design Physical Property Measurement System (PPMS). 
We used the same MPMS to measure the magnetization 
of the Fe-Cr multilayers as a function of external fields 
at temperatures down to 2 K. 

III. RESULTS AND DISCUSSION 

Figure 1 shows the magnetoresistance vs. external field 
H (kOe) for a typical Xe-ion sputtered Fe-Cr multilayer 
sample at 10 and 300 K. The MR becomes constant at 
a saturation field i?sat (vertical arrow) around 13 kOe 
with typical values of 21% at 10 K. These values com- 
pare favourably with 30% and 40% obtained in DC mag- 
netron sputtered and MBE grown samples, respectively. 
Hysteresis in the MR was negligible as the magnetic field 
was swept from (0 -> 20 ^ ^ -20 0) kOe. 



A. Temperature dependence of the electrical 
resistivity 

Figure 2 shows the electrical resistivity (p) vs. temper- 
ature (T) for sample 1 at several values of the external 
magnetic field from to 12 kOe. The saturation field 
(Hsg^t) has a weak temperature dependence, namely, it 
decreases with increasing temperature. This is clear from 
the p{T) curves, say, at 10 and 12 kOe. They are closer 
to each other at higher temperatures. 

To interpret the temperature dependence of the elec- 
trical resistivity, p{T), at an intermediate field between 
H = and H = TJgat is not simple. In this magnetic field 
region the material has both antiferromagnetic and ferro- 
magnetic fractions. Instead, for H > -ffgat, the alignment 
of the spins of each Fe layer is parallel to the direction 
of the external field giving rise to a fully ferromagnetic 
state. It is well known in crystalline bulk 3d metals and 
alloys that the electron-phonon scattering contribution 
to the electrical resistivity from psd dominates over pss 
due to the overlap of the s and d bands at the Fermi level. 
Specifically, for Fe the density of states of the 3d^ ma- 
jority band at the Fermi level is rather large compared 
to those of the s bands. The resistivity, psd, is given by 
the Bloch- Wilson formula, jl^ The "spin-disorder resis- 
tivity" coming from the electron-magnon (spin- wave con- 
tribution) scattering is well described by a relatively 
small term varying as in ferromagnets like Fe, Co, 
and Ni. Putting all these contributions together along 
with the residual resistivity po one can write, assuming 
Mathiessen's rule 

P(r,H = H..., = „ + ^(_) I 

+ BT^ , (5) 

where the second term is the Bloch- Wilson contribution 
Psd and the third term is the small electron-magnon con- 
tribution increasing with temperature due to thermal ex- 
citation of magnons. 

Taking the Debye temperature Od — 420 K for Fe-Cr 
multilayers we have fitted the data for five samples 
at their respective i?sat to Eq. (5) using a 3-parameter 
least-squares fit program which also evaluated the inte- 
gral numerically at each iteration. Excellent fits were ob- 
tained for all the samples with correlation coefficients of 
0.999995 and values of the normalized consistent with 
the number of degrees of freedom and error estimates. 
Figure 3 shows p vs. T data (points) for three samples 
from 5 to 300 K at their respective i?sat- The solid lines 
are the best- fit curves to Eq. (5). It is found that the 
value of the coefficient of the magnetic scattering term, 
B, averaged over all the 5 samples is (4 ± 1) x 10^^ pH. 
cm compared to 1.5 x 10~^ pfl cm K^^ in bulk fer- 
romagnets (Fe, Co, Ni). This higher value of B may be 
related to the fact that the resistivity of these Fe-Cr mul- 
tilayers at 300 K is about five times larger than that of 
bulk iron. 
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The fits of the data of Fig. 3 to Eq. (5) without the 
magnetic {BT^) term are distinctly inferior to those with 
the magnetic term. The values of are typically 6 times 
larger and the correlation coefficients poorer for the fits 
without the magnetic term. The deviation of the actual 
data from the best-fit values (residuals) is plotted in Fig. 
4 as a function of temperature for sample 2 for both the 
fits. The deviation is much less (< 0.1 nD, cm in 40 jiQ. 
cm) and more random for the fits with the magnetic term 
than without it. Addition of a Bloch-Griineissen (BG) 
term (pss), which has a dependence, or replacing the 
Bloch- Wilson term by the BG term makes the fit much 
worse. 



B. Temperature dependence of the 
magnetoresistance 

At H = {] these multilayers are ideally in an antifer- 
romagnetic state where the neighbouring ferromagnetic 
Fe layers are all antiferromagnetically coupled resulting 
in a higher resistivity. Actually there may be pinholes 
through the Cr spacer layer directly coupling the Fe lay- 
ers ferromagnetically instead. Let us define an antiferro- 
magnetic fraction (AFF) as 



AFF(i/) = (1 - M{H)/Ms) X 100% 



(6) 



where Ms is the magnetization measured at iJ = Hg^t 
and M{H) is the magnetization when the field is reduced 
from saturation to H , all at 5 K. Our M{H) measure- 
ments on these samples show that the AFF is typically 
80% Sit H — Q. As H is increased the Fe layers gradually 
turn their magnetization in the direction of the external 
field reducing the AFF and hence the resistivity. Fi- 
nally the AFF reduces to zero (fully ferromagnetic align- 
ment) and the resistivity and hence the GMR saturates 

'dt H — Hsat- 

We define Ap(T) = p{T, H = 0)af- p{T, H = iJ,at)FM 
as the difference in resistivity at a given temperature 
T between the AF {H = 0) and the FM {H = TJ^at) 
states, both assumed ideal. This ^p{T) is primarily due 
to the additional spin-dependent scattering (both bulk 
and interface) in the antiferromagnetic (AF) state. It is 
assumed here that the residual resistivity and the inter- 
band s-d scattering (dominant for 3d metals and alloys) 
do not depend strongly on magnetic fields. Figure 5 plots 
Ap(T) vs. T data (stars) for samples 1 and 3. Thus the 
additional spin-dependent scattering (resistivity) in the 
AF state decreases with increasing temperature. Just 
like the magnetic field aligns the spins in different Fe lay- 
ers reducing the AFF (gradually bringing ferromagnetic 
order in its place) and produces a negative magnetore- 
sistance, here temperature reduces the antiferromagnetic 
order (potentially bringing down the AFF) and hence 
Ap(T). It is seen from Fig. 5 that Ap varies as at 
low temperatures and is roughly linear at higher temper- 
atures. 



Singh et al. Q had worked out the reduction in an- 
tiferromagnetic order due to thermal excitation of spin 
waves in highly anisotropic antiferromagnets with weak 
interlayer coupling between the antiferromagnetic planes. 
This theory has been extended in the present case where 
each Fe layer is ferromagnetic but coupled antiferromag- 
netically to the neighbouring Fe layers due to the RKKY 
interaction. In terms of the planar and interlayer ex- 
change energies Jp and Jz, respectively, the sublattice 
magnetization m{T) at temperature T is given in the 
Appendix (Eq. (13)). Assuming Ap(T) ~ rn{T), we 
have the relation 



Ap(r) 

Ap(0) 



1 T p , 

= 1 2l~ / "9^ In 

Ti" Jp Jo 



1 



(7) 



where m(0) = 1. The expression on the right differs 
only insignificantly from the corresponding expression 
for m{T)/m{0) obtained earlier Q for the anisotropic 
antiferromagnet (cos^ instead of cosq^), where it was 
shown to fall off as at low temperatures (T ^ J^), 
crossing over to an approximately linear (T In T) fall off 
at high temperatures (T ^ J^). We note that this ex- 
pression is relatively unchanged when ferromagnetic do- 
mains are included. 

We numerically evaluated the integral in Eq. (7) and 
used a three parameter least squares fit program to fit 
the data of Fig. 5. The resulting best-fit curves, shown 
by solid lines in Fig. 5, yield values of consistent 
with the experimental errors and a correlation coefficient 
of 0.9999. Estimates of Jp and are found from the 
above fits. They are (230 ± 20) K and (70 ± 20) K, 
respectively. The value for Jp is well below the Curie 
temperature (1040 K) for bulk iron but could be closer 
to the unknown Curie temperature of 20A thick Fe films 
sandwiched between Cr spacer layers. The value for is 
however satisfying, close to the recently identified glass 
temperature, Tg=lAO K, of an antiferromagnetic glassy 
phase that coexists with GMR in similar multilayer films. 
Irreversibilities in this glassy phase have been shown to 
arise from the same interlayer coupling that drives the 
antiparalUel alignments in GMR. 

Consistent with the above model are experimental 
data taken by us (not shown) and others [^| in which 
the saturation fields are studied as a function of spacer 
layer thickness. For three ion-beam sputter deposited 
10-layer samples of Fe(20A)/Cr(c?cr) with different Cr 
spacer layer thickness (dci ) we have observed that as c?cr 
increases from 8 to 12A, the saturation field (i?sat) de- 
creases from 10 to 5 kOe. As Jz decreases with increasing 
dcr, smaller external fields are necessary to break the an- 
tiferromagnetic coupling between the Fe layers. 
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C. Temperature dependence of the 
magnetoresistance in intermediate fields 

Following the work of Aliev et al. M we have fitted 
our data for samples 1,2, and 3 to Eq. (4) with Ps{T) = 
p{T,H) — p(T, iJsat)- The data points along with the 
least-squares fit curves are shown in Fig. 6 for sample 
1 for _ff=0,3,4,6, and 8 kOe. Excellent fits are obtained 
for all the samples with values of consistent with the 
experimental error, correlation coefficients ~ 0.999 
and small errors in the fitting parameters b, c and a. The 
fits are, however, better for smaller fields. Figure 7 shows 
a of Eq. (4) vs. H/Hsat for all the three samples. The 
solid lines are just guides to the eye. The shape of the 
curve is rather similar to the results obtained by Aliev et 
al. 11 (summarized at the end of INTRODUCTION and 
Fig. 3 of Ref. [8]) for MBE grown samples. However, 
we find some differences, like a in our work is typically 2 
for H/Hsat < 1/3, becoming ~ 1 for H / Hsat — 2/3 and 
decreasing at still higher fields. This implies that ps vs. 
T curves (Fig. 6) are quadratic in lower fields and linear 
around H/Hsat — 2/3 instead of ~ 1 as in the work of 
Aliev et al. % 

It is found from Fig. 6 that the intercept h decreases 
with increasing applied fields. This is simply due to the 
fact that the antiferromagnetic fraction (AFF) decreases 
with increasing field. If we plot 6 as a function of our 
measured values of AFF (%) for samples 1 and 3 we find, 
as shown in Fig. 8, that b increases with the AFF in a 
monotonic fashion (both decreasing with H). This is a 
logical conclusion since as iJ ^ 0, the AFF attains its 
maximum value giving the highest resistivity in the ide- 
ally AF ground state. It is to be noted that the decrease 
of h with H and the decrease of Ap with T (Fig. 5) have 
a common origin. It is the decrease of the AFF brought 
about by H and T, respectively. 

D. T vs. i^sat phase diagram for the AF-FM 
transition 

Figure 9 shows the low-field magnetoresistance of sam- 
ple 4 (argon-ion sputtered) vs. external field H for the 
parallel orientation {H in the film plane). It is amply 
clear that the MR ~ in contrast to the findings of 
Aliev et al. § (Fig. 1(a) of Ref. [8]) who found that 
the MR is linear in H for their MBE grown samples. We 
found no linear region in the MR vs. H curve even at 
higher fields until the saturation field (i?sat) of 2-3 kOe 
was reached. As a matter of fact sample 4 (Fig. 9) re- 
fiects an S-shaped curve having points of inflection. How- 
ever, for the perpendicular orientation {H perpendicular 
to the film plane) we find that the MR again goes as 
in agreement with the findings of Aliev et al. |^ 

The observed dependence at low external fields is, 
in fact, expected from simple energy considerations, as 



argued below. The antiferromagnetic ground state of the 
multilayer is characterized by the sublattice magnetiza- 
tion m = (m^ — mB)/2, which takes into account the 
antiparallel orientation of the spin polarization in alter- 
nating Fe layers A and B. The direction of m is arbitrary 
in the ideal isotropic situation. When a small in-plane 
magnetic field is applied, the sublattice magnetization m 
aligns itself perpendicular to the direction of the field. 
This is the lowest-energy configuration as it allows for 
energy gain in all Fe layers due to twisting of spins in 
the field direction. If the twist angle is 9, assumed small, 
then the energy gain is mHsmO ~ mH9. The twist- 
ing also costs energy Jzm^(l — cos 20) w 2Jzm^0^ due 
to loss of antiferromagnetic exchange energy at the layer 
interfaces. Minimizing the net energy change yields the 
optimum twist angle, 9{H) = H/AJm, as proportional 
to the field. Now the reduction in the sublattice mag- 
netization or the antiferromagnetic fraction (AFF), and 
therefore the decrease in resistivity, due to this twist is 
m(l — cos9{H)), which goes as for low fields. 

Figure 10 shows the T vs. i?sat phase diagram for 
sample 4 (argon-ion sputtered, same as that of Fig. 9) 
and sample 1 (xenon-ion sputtered) . Both are in parallel 
orientations. Here 7Jsat(T) is the field at which the MR 
becomes field-independent, i.e., the field-induced AF to 
FM transition is complete. The values of Hsat{T = 0) 
are ~ 3 kOe for sample 4 and ~ 11.5 kOe for sample 1. 
Figure 10 is very similar to Fig. 1(b) of Ref. [8] on MBE 
grown samples for the parallel orientation. 

IV. CONCLUSIONS 

The temperature dependence of the electrical resistiv- 
ity and magnetoresistance has been studied in ion-beam 
sputtered Fe-Cr multilayers. Typical in-plane negative 
giant magnetoresistance (GMR) is 21% at 10 K saturat- 
ing at around 1 tesla. Here each Fe layer is ferromag- 
netic but coupled antiferromagnetically in zero field to 
the neighbouring Fe layers due to the RKKY interaction. 
This gives rise to a high resistance. An external magnetic 
field aligns the spins in different Fe layers producing a 
ferromagnetic alignment beyond Hgat which reduces the 
electrical resistance. The electrical resistivity in the fully 
ferromagnetic state {H = iJgat) between 5 and 300 K 
has been interpreted as the sum of a residual resistivity, 
electron-phonon s-d scattering and spin-disorder resistiv- 
ity ( Eq.(5)). The latter has the same order of magnitude 
as in crystalline Fe. 

We have calculated the sublattice magnetization m{T) 
of these Fe-Cr multilayers in terms of the planar and in- 
terlayer exchange energies (Eq.(13) of APPENDIX). The 
additional spin-dependent scattering in the antiferromag- 
netic state at H = 0, defined by Ap{T) = p{T,H = 
0)af — p{T,H = Hsa,t)FM, is obtained from the experi- 
mental data by assuming that the residual resistivity and 
the electron-phonon scattering are roughly independent 
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of the field. The decrease in Ap{T) with increasing tem- 
perature from 5 to 300 K is explained as arising from 
the reduction in the antiferromagnetic order due to the 
thermal excitation of spin waves, i.e., Ap(T) ~ m{T). 
Mattson et al. [^, on the other hand in dc magnetron 
sputtered Fe-Cr superlattice, found Ap{T) decreasing as 
at temperature below lOOK (Eq. (3)) in contrast to 
our fits (Eq.(7)) over a much wider temperature range. 

From our data at intermediate fields {0 < H < i/gat), 
the spin-dependent part of the electrical resistivity, de- 
fined as p,{T) = p(r, iJ) - p{T,H > H,,t), fits very 
well to Eq. (4). We find that a is typically 2 for 
H/Hsat < 1/3, becoming ~ 1 for H/Hsat — 2/3 and then 
decreasing further at still higher fields. The decrease of 
the intercept b with increasing H and that of Ap{T) with 
increasing T (Eq.(7)) are due to the decrease of the an- 
tiferromagnetic fraction (Eq. (6)) with increasing H and 
T, respectively. Very similar conclusions were reached by 
Aliev et al. |g] in MBE grown Fe-Cr multilayers. 

Finally, we have also obtained the T vs Hgat phase 
diagram for the field-induced AF to FM transition. 



V. APPENDIX 

To obtain the magnon energies in the multilayer sys- 
tem, the following simplified Hubbard model is consid- 
ered on a three-dimensional lattice consisting of a stack 
of layers in the z direction 

+ ^ ■ (8) 



Here the planar band energy ep{kx, ky) together with the 
correlation term describes the ferromagnetic layers, while 
the inter-layer hopping term tz-, which connects sites i to 
nearest-neighbour sites i ± S in the neighbouring layers, 
represents the AF exchange coupling between layers. We 
divide the multilayer system into two sublattices with 
alternating A and B layers, and consider a ground state in 
which the A-layers have spin polarization (nj — n|) = +m 
in the +z direction while B-layers have spin polarization 
(nj — n\) = — m in the — z direction. The sublattice 
magnetization m, a dimensionless quantity, measures the 
AF order parameter in the multilayer system. 

In this two-sublattice basis, and in the HE approxima- 
tion, the Hamiltonian reduces to 



ep(k) - <tA 



ep{k)+aA 



(9) 



where Ok and 6k are the Fourier transforms of the elec- 
tronic annihilation operator Ci , defined on the two sublat- 
tices A and B, respectively. Here 2 A = mil and the sub- 
lattice magnetization m is determined self-consistently. 



For simplicity we consider the strong correlation limit 
in which at the HE level m « 1. The interlayer band 
energy ez(k) = — 2tz cos kz mixes the two ferromag- 
netic bands, and hence the quasiparticle band energies 
i?(k) — ep(k) ± y/ A'^ + ez(k)^ have a mixed character 
with features of both the ferromagnetic and antifer- 
romagnetic ground states, [ p^ As the planar (ferromag- 
netic) band energy ep(k) appears on the diagonal, the 
eigenvectors of the Hamiltonian matrix in Eq. (9) are 
unchanged from the AF case, jl^ 

Evaluation of the magnon propagator x ^(Q"^); i^i" 
volving transverse spin operators {S~ , S^), and repre- 
senting transverse spin fluctuations about the Hartree- 
Fock ordered state, has been described earlier in the ran- 
dom phase approximation (RPA) for both ferromagnetic 
jlSf and antiferromagnetic ground states. For the 
multilayer system the magnon propagator is obtained as 



JpQp - i 



J. 



-Jz COS Qz 



-Jz COS Qz 

Jz + Jpql + 



1 



(10) 



for small planar momentmn qp = (qx,qy). Here Jz — 
2t\l A K, 4:t1/U is the exchange energy characterizing 
the antiferromagnetic coupling between layers, and Jp, 
the magnitude of which depends on details of the planar 
band energy ep(k), plays the role of the planar exchange 
energy. The magnon energy Wq is given by 



Jl COS^ qz 



(11) 



which has the right limiting behaviour yielding the anti- 
ferromagnetic magnon energy Jz\/l — cos^ Qz as Jp ^ 0, 
and the ferromagnetic magnon energy JpQp as ^ 0. A 
completely different starting point in terms of a Heisen- 
berg spin model for the multilayer system, with planar 
and interlayer exchange energies Jp and Jz, would yield 
the same result. 

Going over now to the thermal excitation of magnons 
in the multilayer system, the change 6m{T) = m{T) — 
r7i(0) in the sublattice magnetization at finite tempera- 
ture T is obtained by considering both the advanced and 
retarded modes in the spin-fluctuation propagator with 
appropriate Bose weights. After subtracting out the zero- 
temperature (quantum fluctuation) part the reduction in 
the sublattice magnetization is obtained as 



5m{T) 



2nqpdqp 
(27r)2 



dqz 
27r 



Jz 



JpQp 



- 1 

(12) 



Here the upper limit of integration for the qp integral has 
been taken as oo for convenience, which is valid at tem- 
peratures low compared to Jp, as the high-energy modes 
have exponentially small weight. Integration over the 
planar momentum qp finally yields 



6 



m{T)=m{0)~ r dq. In 

Jp Jo 



-^(l-COs2 g,)l/2 

(13) 



FIG. 1. Magnetoresistance vs. external field H (kOe) ori- 
ented parallel to the layers for a Xe-ion beam sputtered Fe-Cr 

multilayer sample (sample 2) at 10 and 300 K. The MR sat- 
urates around 13 kOe (//gat) and has a typical GMR of 21% 
at 10 K. 
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FIG. 2. Electrical resistivity (p) vs. temperature (T) for 
sample 1 at several fields between and 12 kOe. The curves 
are closer to each other at higher temperatures indicating that 
the saturation field Hsat decreases with increasing tempera- 
ture. 

FIG. 3. Electrical resistivity (p) vs. temperature (T) data 
(points) from 5 to 300 K for samples 1, 2 and 3 at their 
respective Hgat . The solid lines are the excellent least-squares 
fitted curves for fits to Eq. (5) which includes lattice and 
magnetic scattering contributions. 
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FIG. 4. The deviation of the actual data from the best fit 
values is plotted as a function of temperature (T) for sample 2 
for fits with and without the magnetic term. The deviation 
is much less (< 0.1//f2 cm in 40 /xf2cm) and more random for 
the fits with the magnetic term than that without it. 

FIG. 5. Ap(T) = p{T,H = 0)af - p{T,H = H,^^,fm 
vs. temperature(T) data (stars) for samples 1 and 3. This 

additional spin-dependent resistivity in the AF state decreases 
with temperature. Ap varies as at low temperatures and 
roughly linearly at higher temperatures. The solid lines are 
the least-squares fitted curves for fits to Eq. (7). 

FIG. 6. ps{T) = p{T,H) - p(r,_ffsat) vs. temperature (T) 
data (points) for sample 1 at H=0, 3, 4, 6 and 8 kOe. The 
solid lines are the least-squres fitted curves for fits to Eq. (4). 

FIG. 7. a of Eq. (4) vs. H/Hsat for all the three samples. 
The solid lines are just guides to the eye. 

FIG. 8. b of Eq. (4) vs. AFF(%) for samples 1 and 3. 6 is 
found to increase monotonically with the AFF (both increas- 
ing with H). 



FIG. 9. Magnetoresistance vs. external field H(kOe) at low 

fields for sample 4(Ar-ion sputtered) for the parallel orienta- 
tion (H in the film plane). The temperatures are between 5 
K and 205 K (every 20 K) and 300 K. Clearly the MR ~ H"^ 
at lower fields. 
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FIG. 10. T vs. Hsax phase diagram for samples 4 
(Ar-ion sputtered) and 1 (Xe-ion sputtered). Here ffsat is 
the field at which the MR becomes field independent, i.e., the 
field-induced AF to FM transition is complete. The values of 
Hs3x(T = 0) are 3 kOe for sample 4 and 11.5 kOe for sample 
1. 
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